Sr2Cu(P04)2: A real material realization of the ID nearest neighbor Heisenberg chain 
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We present evidence that crystalhne Sr2Cu(P04)2 is a nearly perfect one-dimensional (ID) spin- 
1/2 anti- ferromagnetic Heisenberg model (AHM) chain compound with nearest neighbor only ex- 
change. We undertake a broad theoretical study of the magnetic properties of this compound using 
first principles (LDA, LDA-I-U calculations), exact diagonalization and Bethe-ansatz methodologies 
to decompose the individual magnetic contributions, quantify their effect, and fit to experimental 
data. We calculate that the conditions of one- dimensionality and short-ranged magnetic interac- 
tions are sufficiently fulfilled that Bethe's analytical solution should be applicable, opening up the 
possibility to explore effects beyond the infinite chain limit of the AHM Hamiltonian. We begin 
such an exploration by examining some extrinsic efi'ects such as impurities and defects. 
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I. INTRODUCTION 



Interest in magnetically low-dimensional systems be- 
gan with the advent of quantum mechanics and the devel- 
opment of spin-spin interaction models to explain mag- 
netic behavior. Though deceptively simple, early mod- 
els such as Isingi or Heisenberg^ cannot be solved for 
the general case, requiring either low spin or low spa- 
tial dimensionality to obtain analytical solutions. Ex- 
act solutions for some specific cases, such as Onsager's 
solution^ to the two-dimensional (2D) Ising model or 
Bethe's solution'* to the spatially one-dimensional (ID) 
s=l/2 Heisenberg chain, inspired general interest in the 
theoretical properties of magnetism in low dimensions, 
and resulted in a variety of predictions for exotic phys- 
ical behaviors. Some of the most interesting proper- 
ties to arise from the study of restricted-dimension mag- 
nets are due to the dramatic effect of quantum spin 
fluctuations^. These are intimately involved in the emer- 
gence of novel ground states and excitation spectra, 
and, in recent decades, have become a leading candi- 
date as the pairing mechanism for electrons in quasi-2D 
high temperature superconductor^, thus providing the 
field with a more practical aspect. Mermin and Wag- 
ner demonstrated'' that strong spin fluctuations suppress 
magnetic long range order (LRO) until T=0 in isotropic 
spin systems in ID and 2D, regardless of the interaction 
strength (exchange) between neighboring spins. Since 
real physical compounds are 3D by nature, the continu- 
ing effort to experimentally verify such predictions and 
probe the nature of ID magnets is aimed at finding mate- 
rials where magnetic interactions proceed predominantly 
along one direction. A measureable gauge of success or 
failure along these lines can be obtained through the ra- 
tio fcsT/v/ Ji, which compares ordering temperature of a 
Neel state to the magnitude of exchange between neigh- 
boring spins. A perfectly ID system would give ksT^ / Ji 
— 0. Additionally, experimentally observed behaviors 
should conform to theoretical predictions where the con- 
ditions of the model are met. 



Recently, very low temperature measurements^ on the 
spin-1/2 compound Sr2Cu(P04)2 identified the onset of 
magnetic LRO at Tjv= 0.085 K, which, in combina- 
tion with the extracted exchange constant, yields a ra- 
tio ksTpf/Ji = 6x10""*. This is nearly a full order of 
magnitude less than the ratic^ for prototype ID magnet 
Sr2Cu03, which has kBT^/Ji = 2.5x10"^. We can esti- 
mate the remaining interchain coupling, J', by adopting 
an effective 3D chain model with z±—4: nearest neighbor 
chains as in Ref. ITol : 
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where kAPM is the magnitude of the AFM wavevec- 
tor. All signs, therefore point to an extremely high de- 
gree of one-dimensionality that should justify the use 
of Bethe's exact analytical solution to the ID spin-1/2 
AHM in a wide temperature range J' -C fcsT, pro- 
vided that indeed only nearest-neighbor interactions are 
present. Previous studies^^ have shown that this con- 
dition is not satisfactorily fulfilled by the leading ID 
spin-1/2 AHM candidate, Sr2Cu03. For completeness, 
another candidate for a ID- AHM system should be men- 
tioned: the linear charge transfer salt [3,3'-dimethyl-2,2'- 
thiazolinocyanine]-TCNQi2ii^. 

Here, we undertake an extensive theoretical analy- 
sis of Sr2Cu(P04)2, employing first principles density 
functional theory calculations, an exact diagonalization 
scheme, and finally, a Bethe-Ansatz derived expression 
for the magnetic susceptibility fit to experimental data. 
Our results are in remarkable agreement with one an- 
other and with experimental observations. We show 
that Sr2Cu(P04)2 is indeed extremely ID and that fur- 
thermore, the second-neighbor interactions are vanish- 
ingly small, eliminating any complications due to in- 
chain frustration. We claim, therefore, that this com- 
pound is the best realization of a nearest-neighbor only 
Heisenberg spin 1/2 chain known to date and will likely 
be a valuable tool for extracting intrinsic effects be- 
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FIG. 1: (Color online) Upper panel The structure of 
Sr2Cu(P04)2, showing the isolated CUO4 plaquettes as 
squares with O ions located at the corners, and the PO4 units 
as tetrahedra with P located in the center. Lower panel A top 
view of the spin chain plane with the various hopping paths 
labelled. Hopping to the nearest out-of-plane neighbor, tx, is 
not shown. 



yond the Bethe-ansatz i.e. effects not contained in the 
Heisenberg Hamiltonian, such as Dzyaloshinsky-Moriya 
interactionsi^, ring exchange processes, or spin-lattice 
coupling effects. Additionally, extrinisic effects due to 
sample imperfections such as defects, magnetic impuri- 
ties or the presence of alternate phases can be quantita- 
tively explored, and we make a preliminary investigation 
of some of them. 



II. FIRST PRINCIPLES AND TIGHT-BINDING 

First principles density functional theory (DFT) band- 
structure calculations were performed using a full- 
potential local orbital code, FPLOi^, with the following 
valence states: Sr (4s,4p,4d,5s,5p), Cu (3s,3p,3d,4s,4p), 
P (2s,2p,3s,3p,3d), O (2s,2p,3d). The structure, lattice 
constants, and atomic positions (see Fig. ^ were taken 
from experiment ^^-iii, a = 11.515 A, 6 = 5.075 A, c — 
6.5748 A. Sr2Cu(P04)2 and isostructural Ba2Cu(P04)2, 
which we calculate as a comparison material, are com- 
posed of isolated CUO4 plaquettes surrounded by PO4 
tetrahedra. The plaquettes are spaced evenly along the 
b axis, forming chains that are staggered with respect 
to one another in the a crystal direction. The planar 



CUO4 units are tilted with respect to the a — b crystal 
plane. Based on the geometry of these two systems, we 
enumerate five specific interactions between plaquettes, 
and therefore between spins localized to these plaquettes, 
that may be necessary to describe the electronic and mag- 
netic structures. These interactions are labelled in Fig. 
n as various hopping parameters to be later included in 
a tight-binding model. 

The paramagnetic band structure of Sr2Cu(P04)2 
(Fig. 12) shows a single, isolated, half-filled band, derived 
predominantly from the Cu 3dx2_y2-0 2p^ molecular pla- 
quette orbital, crossing the Fermi energy. In reality, the 
system is antiferromagnetic (AFM) and insulating, but 
we will follow the standard procedure of importing hop- 
ping parameters from the " uncorrelated" paramagnetic 
system to the Hubbard model which then maps onto the 
Heisenberg model. Magnetism and correlation effects can 
also be added at the DFT level with, as we will show, 
very consistent results. The ID character of the sys- 
tem is qualitatively obvious from the nearly dispersion- 
less bands in directions perpendicular to the magnetic 
chains (S-Y and F-Z) and from the characteristic loga- 
rithmic divergences in the density of states (DOS) near 
the band edges. To quantitatively compare microscopic 
magnetic interactions, we fit a tight-binding (TB) model 
to our paramagnetic band structure and calculated the 
individual exchange constants between various CUO4 pla- 
quette spins using Jij = 4t| /Ue// with Ue//=4.5 eV. 
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FIG. 2: (color online) The bandstructure of Sr2Cu(P04)2, 
showing the single metallic band well separated from all oth- 
ers. The X — S and Y — Y directions are along the magnetic 
chain, while the S — Y and V — Z directions are perpendicular 
to the chain. The inset shows the total and orbitally resolved 
DOS for the single band. In the lower panel, a blow-up of this 
band is shown with tight-binding eigenvalues superimposed to 
show the remarkable reproduction of the dispersion. 
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We find this to be a reasonable choice for Ue/ / because 
of poor screening in this geometry and because of the 
smah inter-plaquette repulsion. The hopping parame- 
ters included in the model are shown schematically in 
Fig. ^ with the numerical values and derived superex- 
change constants given in Table H The resulting TB dis- 
persion, which uses only 5 fitting parameters, is indistin- 
guishable from the full-potential calculation (Fig. [flower 
panel), indicating that further interactions can be safely 
ignored. The ratio of the strongest in-chain coupling to 
the strongest interchain coupling is Ji/J!'^ ^ 70 and the 
ratio of first to second neighbor in-chain coupling is J1/J2 
~ 700 for Sr2Cu(P04)2. Identical calculations based on 
the band structure of Ba2Cu(P04)2 (not shown) yield 
similar results with somewhat more inter-chain coupling 
but less second neighbor in-chain coupling. Both systems 
can therefore be considered as strongly one-dimensional, 
with Ba2Cu(P04)2 slightly less so. Naturally, the choice 
of Ue/ / is simply a best estimation and results for J will 
vary slightly with this choice, while the ratios will remain 
constant. 

TABLE I: Tight-binding hopping parameters (in units of 
meV) and derived exchange constants (in units of K) for 
A2Cu(P04)2, A=Sr,Ba. See Fig. E]for the relationship of 
the hopping parameters to the structure. 
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0.008 


1.09 


0.228 


0.03 



2p contributions in addition to Cu 3d) increasingly lo- 
calizes the Cu 3c? electrons. For a range of Ud between 
6 eV and 9 eV, we find that Jeff varies from 261 K to 
160 K. Since is a local quantity and since the Cu-0 
bond distance in Sr2Cu(P04)2 is only 1% different than 
in the plaquettes of the widely studied high Tc precur- 
sor systems, we adopt the commonly used value of Urf = 
8 eV. This corresponds to a value of 190 K for Jeff - Note 
that since the LSDA (and LSDA-I-U) energy differences 
include contributions from all exchange processes in the 
system, Jef / cannot in general be considered as either 
purely superexchange or purely ID. However, comparison 
with the individual superexchange parameters derived 
from the TB fit shows that both assumptions, in this 
case, are quite valid. The nearest neighbor in-chain TB 
exchange constant has a value Ji = 187 K, in exceptional 
agreement with the Je// value of 190 K, indicating that 
the next-nearest neighbor interactions, FM exchange pro- 
cesses and residual 2D and 3D interactions must therefore 
be extremely small. Of course both values can be made 
to vary somewhat by choosing Ud and Ue// differently, 
thus affecting the agreement as well. We expect that our 
calculated value of J will be larger than the experimental 
value, as it is well known that the band dispersion from 
which we derive t and subsequently J is generally exag- 
gerated by the LDA. Indeed, the experimentally derived 
value of J is 143 K, in good but not perfect agreement 
with our calculations. We emphasize that any renormal- 
ization of the hopping parameters stemming from effects 
outside the LDA will cancel in the ratio {J1/J2) so that 
the precise calculated value of J has, in any case, no bear- 
ing on the establishment of the compound's pronounced 
magnetic one-dimensionality and short-ranged magnetic 
interactions that are the primary aim of our first princi- 
ples study. 



The energy difference between FM and AFM ordered 
spin configurations can be calculated using the local spin 
density approximation (LSDA)i& which allows for sepa- 
rate spin-up and spin-down densities. Since the LSDA 
is known to badly underestimate the onsite Coulomb in- 
teraction in localized systems, we applied the LSDA-fU 
methodology to better account for the correlated Cu 3d- 
orbitals, using the fully localized limit schemei^ to cor- 
rect for double-counting terms. We map a classical Neel 
state and a classical ferromagnetic state onto the Heisen- 
berg spin model, including only ID nearest neighbor in- 
teractions. Comparing the resulting model energy dif- 
ference to the LDA energy difference between FM and 
AFM states (per spin), we derive an effective exchange 
constant, Je//, in the following way: 

H — ^ JijSi-Sj-, Epm—Eafm = 2Je//|sp, s — 1/2, 

(2) 

As expected, the energy difference, and therefore Je//, 
decreases as U^ (not to be confused with the consider- 
ably smaller one-band parameter Ue// that contains O 



III. EXACT DIAGONALIZATION 

We perform an exact diagonalization calculation us- 
ing ten sites along two staggered AFM magnetic chains 
(20 sites total) for calculating thermodynamic properties 
such as specific heat and magnetic susceptibility, and us- 
ing 36 total sites for obtaining the ground state proper- 
ties. We compare calculations including the three largest 
exchange interactions, Ji, J2 and J?e, listed in Table 
^ to calculations using only Ji . Since in our model each 
chain has only one neighboring chain, we also perform 
a calculation in which Jl^ is doubled to account for the 
existence, in reality, of two neighboring chains. In the 
case of specific heat, there is no discernable difference 
between any of the three curves using these different pa- 
rameter sets. For the calculation of x, the curves are 
identical for the majority of the temperature range ex- 
plored (0 < khT/Ji < 3), but a barely visible difference 
occurs near the peak of the curve (see Fig. |3Jl. The max- 
imal difference occurring between two curves is 0.9%, at 
about ksT/ Ji=0.65. The ground state calculations of 
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FIG. 3: A 20 site exact diagonalization calculation of the 
susceptibility per spin of Sr2Cu(P04)2. The inset shows a 
blow-up of the region where the curves have the greatest dis- 
crepancy - 0.9%. Curves for the ground state correlation of 
spins and for the specific heat per spin show even less devia- 
tion. 
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FIG. 4: Bethe-Ansatz fits to experimental data with correc- 
tions for broken chains, impurity phases, and domain bound- 
aries. Note the consistency of the exchange parameter vs. the 
wide spread in values for the chain length. 



the spin-spin correlation, (Si ■ Sj) are again completely 
indistinguishable. Obviously, more distant and weaker 
interactions, such as Jf^ and J± will have even less of 
an effect. From this we conclude that Sr2Cu(P04)2 is 
essentially free of magnetic interactions beyond the first 
neighbor and is an ideal candidate for Bethe-ansatz cal- 
culations, which we now discuss. 



IV. COMPARISON WITH EXPERIMENT 

The basic theory for non-periodic, open AHM chains 
was worked out by Furusaki and Hikihara^" (FH) and 
also by Zvyagin and Makarova^i (ZM). FH considered 
a half-infinite chain with one free chain end applying 
bosonization theory whereas ZM considered finite even- 
membered chains on the basis of a rigorous theory based 
on the Bethe ansatz. Both approaches result, at low 
temperature, in a chain length (L = Na) dependent di- 
verging contribution to the total magnetic susceptibil- 
ity X cx: l/NT\n{ToJi) and to the linear coefficient in 
the specific heat 7 = Cp/T oc l/iVTln'^(To Ji/T), with 
(3=2,4 in the FH and ZM theory, respectively. Since a 
real chain has two ends we multiplied the FH expres- 
sions by a factor of 2. ZM calculated further logarithmic 
corrections which we adopted here to be valid for the FH 
case, too. Then within both approaches we arrive finally 
at the same expression for the chain end contributions 
to X (up to a factor of 3/4). For the logarithmic con- 
stant Tq we used the same value 5.696 as proposed by 
Johnston et al. in the fit expression (fit2) for the bulk 
susceptibility^^, similarly to 5.8 used in Eq. (1). In the 
shown and described fits we have adopted the formalism 
of FH, modified as described above, for both C/T and x- 



TABLE II: The collected results of fits to different sets of 
X(T) data for Sr2Cu(P04)4 and Ba2Cu(P04)4. 



Sr2Cu(P04)2 


J(K) 


N 


g 


Belik^'^ 


144.9 


114 


2.713 


Belik® 


143.9 


337.7 


2.154 


Nath^^ 


152.6 


5 


1.952 




Ba2Cu(P04)2 


J(K) 


N 


g 


Belik^** 


131.8 


48 


2.073 


Nath2^ 


140.1 


5 


2.041 



We fit to both susceptibility and specific heat to data 
taken on the same sample^^, additionally using data from 
a later samplc^"^ for which only x data was available. 
The fit quality for specific heat and susceptibility are 
very similar; results for the latter are shown in Fig. 0] 
We get very good fits throughout the entire tempera- 
ture range and find a reasonably consistent value for the 
exchange parameter, J, despite fitting to samples of dif- 
ferent qualities and to two different measurements (C/T 
and x). It is worth noting that the exchange parame- 
ter derived from fitting to C/T using the data of Ref. 
IT^ yields an exchange parameter of J=134.4 K, which 
is somewhat less than the value derived from fitting to 
x{T) using the same sample. This is likely caused by 
non-magnetic impurity contributions (not accounted for 
in our model) that affect the specific heat but not the 
susceptibility at low temperatures. Since the magnetic 
component of C is oc 1/J, ignoring the non-magnetic con- 
tribution overestimates this term and therefore underes- 
timates the exchange. Unlike the exchange and g values, 
the chain length parameter, iV, varies widely between 
samples. While this is expected for samples of different 
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quality, the variation is surprisinly high and, more impor- 
tantly, the chain lengths resulting from our fits are far too 
small to justify the use of the open chain theories that 
we have employed at very low temperatures. With chain 
lengths of this order, the low temperature region will be 
completely dominated by broken chain physics that re- 
quires different, and as-of-yet undeveloped, formalism. 
Having fit throughout a large temperature range, includ- 
ing regions where broken chain physics is inoperative, we 
feel that the extracted exchange constant, J « 145 K, 
is nonetheless relevant - a belief that is supported by its 
consistency between fits and its similarity to the experi- 
mentally measured value. 



V. DISCUSSION 

The application of a variety of theoretical techniques 
to the problem of magnetism in Sr2Cu(P04)2 convinc- 
ingly demonstrates that the ideal compound is highly 
one-dimensional. The high degree of one-dimensionality 
can be traced back to its unusual isolated CUO4 plaquette 
geometry. Instead of edge- or corner- shared plaquettes 
such as are common in other quasi- ID compounds^iiS^, 
each Cu ion in Sr2Cu(P04)2 is surrounded by four O ions 
not shared by any other Cu ion. This construction virtu- 
ally eliminates the second neighbor in-chain coupling that 
prevents edge-shared compounds such as Li2Cu02 from 
being described via a simple nearest neighbor Heisenberg 
modelSSiSl. Corner shared cuprates such as Sr2Cu03 
have far smaller second neighbor interactions^, of the or- 
der J1/J2 ^ 15, and yet, these must be taken into account 
to get good agreement between model calculations and 
experiment The structure of Sr2Cu(P04)2 along the 
chain is that of an edge shared chain compound with ev- 
ery other unit missing. Conceptualized in this way, one 
can make a correspondence between exchange constants 
in a edge-shared (es) system and those in the isolated 
square plaquette (sp) geometry: J^** Jj^, and Jf'* — > 
Jj^. Since Jf" is known to be vanishingly small in the 
edge-shared geometry, it is clear that the second neigh- 
bor interactions in the square plaquette geometry can be 
expected to be negligible. This may provide some di- 
rectional guidance in the search for new one-dimensional 
compounds: the isolated plaquette arrangement appears 
to be superior to the more common edge- or corner- 
shared structures such that synthesis of new compounds 
with this geometry may prove to be profitable. The tilt- 
ing of the out-of-chain plaquettes with respect to one 
another further suppresses the inter-chain coupling. The 
staggering of plaquettes in neighboring chains slightly in- 
creases the distance between spins, but more importantly, 
gives rise to frustration. As each chain is antiferromag- 
netically aligned by the (relatively) strong first neighbor 
coupling, a given spin finds itself surrounded by four in- 
terchain neighbors, two aligned in one direction and two 
in the other. These staggered, frustrated chains are more 
decoupled from one another than they would be in an- 



other arrangement, e.g. a ladder configuration. 

Provided that Sr2Cu(P04)2 is stoichiometric and 
largely defect-free, it is clear that this compound repre- 
sents the most ID AHM chain so far investigated. These 
conditions are, unfortunately, not reasonably fulfilled by 
current samples. It is interesting in this context to con- 
sider the mechanism by which the compound eventually 
achieves LRO (at = 0.085 K): is it truly the result of 
residual third dimension interactions? Significantly, the 
phenomenologically estimated averaged interchain inter- 
action from Eq. ^is of the same order as the calculated 
Jj_, rather than J^^. Fluctuation induced "order by dis- 
order" coupling could be responsible for the strong re- 
duction of two orders of magnitude within the frustrated 
plane. On the other hand, the interchain couplings are in 
general phenomenally small as calculated by DFT meth- 
ods and even so are likely exaggerated. One alternative 
explanation is that in a system with many broken chains, 
there will be some number of chains containing an odd 
number of spins, with each such chain carrying one un- 
compensated spin- 1/2 electron. The relationship of the 
uncompensated spins to one another is not defined by 
any of our methods and a long range ordering of these 
is not out of the question. It would be interesting to 
see if the ordering temperature remains constant with 
sample quality. Another point to be addressed in the fu- 
ture is the issue of spin-lattice coupling. The Heisenberg 
model itself assumes perfect isotropy in spin-space and 
we have not included any relativistic (spin-orbit) inter- 
actions in our first principles calculations. The neglect of 
these is seemingly justified by the extremely small field 
(H= 4mT) at which the spin-flop transition occurs^, but 
the smallness of this field itself is unusual and a cause 
for further investigation. All of these facts point to the 
high desirability of better samples that can be used to 
disentangle true "dimensionality" effects from behaviors 
due to crystal imperfections. Although the investigation 
of true Heisenberg physics is currently limited by sample 
quality issues, Sr2Cu(P04)2 is theoretically, and poten- 
tially experimentally, the best example of a magnetically 
ID crystal yet studied. 



VI. CONCLUSION 

We have shown that the isolated CUO4 plaquette ge- 
ometry of Sr2Cu(P04)2 gives rise to a nearly perfect ID 
spin-1/2 nearest neighbor only system. We find a ratio 
ksTN / Ji = 6xl0~^, in good agreement with experimen- 
tal finding and show that secondary interactions (2D, 3D 
and next-nearest neighbor) are negligible in terms of cal- 
culated thermodynamic properties. Using the Bethe- 
ansatz solution to the Heisenberg Hamiltonian along with 
additional terms to correct for extrinsic non-crystalline 
effects, we fit the data over a large temperature range 
and derive an exchange parameter of 145 ± 5 K that is 
consistent between samples and between fitting choices. 
We find that sample quality, particularly the existence of 
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numerous broken chains, currently prohibits experimen- 
tal observation of true spin-1/2 AHM physics. However, 
Sr2Cu(P04)2 is truly intrinsically perfectly ID with only 
one exchange parameter, and as better and better meth- 
ods of generating the compound emerge, effects beyond 
Bethe-ansatz can be probed experimentally. 
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